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_  T”*^work’  Wfr}ie  Special  Functions  -.aid  Theln  ApororWti-'ns  *  & 
class  of  rational  approximations  for  the  Generalised  ^ 

p'  ‘  ‘  (Jm(z/  -'-n  be  expressed  in  terms  of  a" /nv, .  or  a  7  ?- 

’  ®orrc spondinj  to  each  form  nrd  a  choice  of  certain  ' 

fcnere  rs  ji_ratienal  aoproxiiaatj  on  for  dr  ( ■?  V  t  n  n 

r-.,of  T  ri  iorjxv(z)  .  j.  c.  p.  Miller  has  ehcvn 

u“~t  ,-s+v(”)a  ra  a  positive  into-'^r  or' -em  1  „  , 

-y.  .......  -  .  ,  •v.hsjcd  or  can  lie  aporoximafs-.i  rv  ■]« 

°-  recursion  torrr.ula  for  /f  {T)<  - 

y-t  4.V.J  1  ,  vn+v^/  ^  pX.Lcd  in  tr.G  b<xck'.rc?.i'J  <-•>■> 

If  tn.s  scheme  is  used  togetfc^OTi  each  of  two  certain  nor-U-ation  "  ’ 
relations,  then  rational  approximations  forffjfj'.  e^r-e  ard w,  , 
appr - ximati ens  arc  identical  those  notrd4L . 


CC..0  it  is  Ito,  that  the  m*wt  error  is  i^ii;nm0ant.  ‘  ‘ 

— %  . ...  ^ 

.  l  ,  i.  a  .)  ... ;  -i  V 


ABSTRACT 


In  my  work,  "The  Special  Functions  and  Their  Approximations,"  a 
class  of  rational  approximations  for  the  generalized  hypergeometric  functions 
was  developed.  Now  Iv(z)  can  be  expressed  in  terms  of  a  qFj_  or  a  ^F^  . 
Thus,  corresponding  to  each  form  and  a  choice  of  certain  free  parameters 
there  is  a  rational  approximation  for  I^(z)  .  J.  C.  P.  Miller  has  shown 
that  Inn-V(z)>  m  a  positive  integer  or  zero,  can  be  approximated  by  use 
of  the  recursion  formula  for  z )  applied  in  the  backward  direction. 

If  this  scheme  iz  used  together  with  each  of  two  certain  normalization 
relations,  then  rational  approximations  for  I^(z)  emerge  and  these  rational 
approximations  are  identical  with  those  noted  above.  The  analysis  leads  to 
&  new  interpretation  of  the  backward  recursion  scheme.  We  also  study  a 
third  case  for  the  evaluation  of  Im+S)(z)  ,  m  a  positive  integer,  by  the 
backward  recursion  process  which  presumes  that  Iv(z)  is  known.  In  each 
instance  a  closed  form  expression  for  the  truncation  error  is  developed 
which  leads  to  a  very  effective  a  priori  estimate  of  the  error.  For  each 
case  it  is  shown  that  the  round-off  error  is  insignificant. 


INTRODUCTION 


0 


4» 


In  my  treatise  on  the  special  functions  [ii  ,  a  class  of  rational 
approximations  for  the  generalized  hypergeometric  function  pFq  was 
developed.  These  approximations  depend  on  a  number  of  free  parameters. 

Since  Iv(z)  can  be  expressed  in  terms  of  a  QF^  or  a  ^F^  ,  there  is  a 
particular  rational  approximation  corresponding  to  each  of  these  hyper¬ 
geometric  forms  and  a  choice  of  the  aforementioned  free  parameters. 

The  idea  of  using  the  recursion  formula  for  I^z)  in  the  back¬ 
ward  direction  to  generate  values  of  Iv(z)  is  due  to  J.  C.  P.  Miller  [2]  . 

It  is  a  very  powerful  tool  and  the  notion  has  created  considerable  interest; 
see  [l,  Vol.  2,  pp.  159-166],  [3,4]  and  the  references  quoted  in  these 
sources.  The  Miller  scheme  together  with  two  certain  normalization  relations 
also  gives  rise  to  rational  approximations. 

In  a  conversation  Jerry  L.  Fields  conjectured  that  the  specific 
rational  approximations  noted  in  the  first  paragraph  are  identical  to  the 
certain  rational  approximations  which  emerge  by  use  of  the  backward  recur¬ 
rence  scheme  noted  in  the  second  paragraph.  In  the  present  paper,  we  verify 
this  conjecture.  In  addition,  we  develop  a  new  interpretation  of  the  Miller 
method.  We  also  study  a  third  normalization  technique  which  is  sometimes 
used  with  the  backward  recursion  scheme.  A  closed  form  analytical  expression 
of  the  error  for  each  case  is  derived.  These  equations  are  valuable  as  they 
lead  to  simple  asymptotic  estimates  of  the  error  which  are  very  realistic 
and  easy  to  apply  in  practice.  It  is  demonstrated  that  the  round-off  error 
is  insignificant.  The  paper  closes  with  some  numerical  examples. 

In  the  main  body  of  the  paper,  we  find  it  convenient  to  deal  with 
the  modified  Bessel  function  Iv(z)  .  The  results  are  valid  for  all  z  in 
the  cut  complex  z-plane  -  rr  <  arg  z  <  tj  and  in  the  cut  complex  v-plane, 

|arg  v  |  <  tt  .  In  this  connection,  we  should  note  that  I_v(z)  =  Iv(z)  if 
v  is  an  integer  or  zero.  Thus  we  suppose  throughout  that  v  is  not  a 
negative  integer.  Actually,  it  is  sufficient  to  have  0  s  arg  z  s  tt/2  in 
view  of  the  definition  of  Iv(z)  .  Also  it  is  sufficient  to  have  R(v)  >  -1  , 
for  if  Ij__v(z)  and  I_v(z)  are  known,  computations  of  I_m_v(z)  , 
m  -  a  positive  integer,  can  be  done  by  use  of  the  recursion  formula  for 
I-m-v(z)  •  All  of  this  not  withstanding,  it  is  convenient  to  restate  some 
of  the  key  equations  to  facilitate  application  of  our  results  to  the  Bessel 
function  Jv(z)  .  This  is  done  near  the  end  of  the  paper. 
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whence 


hji(z)  =  (-)nn!(n+v<-2-6)nXn[l+0(n-1)]  , 


(9) 


4. 


°^)  M*-1)]  >  R^)  >  -1  «  (10) 

and  so,  for  z  and  v  fixed,  R(v)  >  -1,  the  approximation  process  is 
convergent . 


Proof:  Equations  (2) -(5)  follow  from  [l,  Vol.  2,  p.  96]  with 


a=0,f=g=0,p=0,q=l,p1=vH, 
a  =  1  -  6  ,  3  =  v  ,  X  =  v+2-6  ,6=Oor6=l,v=z 

p 

and  z  replaced  by  z  /4  .  Notice  that  the  3F0  series  for  hn(z)  in 
(5)  turned  around  is  the  alternative  form  for  h^z)  in  (6). 

Equation  (8)  follows  from  [l,  Vol.  2,  p.  103]  while  (9)  follows 
from  (6),  see  also  [l,  Vol.  1,  pp.  259-26l],  and  (10)  is  now  obvious. 


Remark:  In  the  proof  developed  in  the  cited  source,  it  was  necessary  to 
suppose  that  R(v)  >  -1  .  Later,  we  present  a  new  formulation  of  the  error 
which  shows  that  v  is  unrestricted  save  that  v  is  not  a  negative  integer. 
So  throughout  this  work  \>  is  arbitrary  except  as  just  indicated.  Compu¬ 
tation  wise,  the  exception  is  no  burden  since  I_n(z)  =  In(z)  . 


Theorem  2.  Both  ijrn(z)  and  hn( z )  satisfy  the  same  recurrence  formula 


Mz)  +  (Ci+XD1)hn.1(z)  +  (C2+XD2)'hn_2(z)  +  C3hn_3(z)  =  0  , 
=  (2n+\-g)(n-X+l)  _  n(2n+\-2)(2n+\-l) 

1  (2n+x-4)(n+X-l)  ’  "  n+X-1  ’ 


_  -(n-l)(n+2X-4)(2n+> -1)  _  -(n-l)(2n+x-2)(2n+X-l)(n+\-3) 

2  (n+x-l)(n+X-2)(2n+X-5)  *  D2  “  (n+X-l)(n+X-2) 
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(11) 


whence 


V*) 


(2/2)V(2Z)n(v+3/2)n0(nu>) 

r(sH-l)(n+2\H-2)nn: 


(19) 


and  for  z  and  v  fixed,  R(\>)  >  -1  ,  the  approximation  process  is  con¬ 
vergent  . 


Proof:  Equations  (13) -(15)  follow  from  [l,  Vol.  2,  p.  96]  with 

a=0,  f=g=0,p=q=l,a1=  ,  Pl  =  2v+l  , 

or=l,0  =  2v,X  =  2\tf-2  ,  y  =  z 

and  z  replaced  by  -2z  .  Equation  (16)  is  equation  (15)  turned  around. 
Equation  (17)  comes  from  [l,  Vol.  2,  p.  103]  while  (18)  comes  from  (16), 
see  also  [l,  Vol.  1,  pp.  133,  259-261].  Thus  (19)  is  at  hand. 


Remark:  See  the  remark  after  Theorem  1. 


Theorem  4.  Both  Sn(z)  and  gn(z)  satisfy  the  same  recurrence  formula 


gn  (z )  +  (E1+2Fi/z)gn_i(z)  +  (E2+2F2/z)gn_2(z)  +  E-jg^z)  =  0 


-(n+v)(n+2v-l) 


1  (n+v-l)(n+2v^-l)  *  F1  n+2v+l 


4n(n+v) 


r  _  ^(p-1)  v  =  -4(n-l)(n+v)(n+2v-l) 

J2  (n+2v)(n+2v+l)  ’  2  (n+2v) (n+2v+l) 


E,  =  _ (n-l)(n-2)(n+y) 

5  (n+2v) (n+2v+l) (n+v-l) 


(20) 


where  n  s  3  . 

Proof:  See  [l,  Vol.  2,  Ch.  12]. 
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BACKWARD  RECURRENCE  SCHEMATA  FOR  GENERATING  Iv(z) 


The  technique  for  generating  I^(z)  by  use  of  the  recurrence 
formula  for  Iv(z)  employed  in  the  backward  direction  is  as  follows.  The 
recurrence  formula 


<Pm 


,  s  2(m+v+l) 

,v(z)  =  “  L  cPmt-l,'j\z)  +  tpm+-2,v(z) 


(21) 


is  satisfied  by 


i{m+u)n 

Vv(z)  and  e  W2) 


Kin+\)(z)  sin(m+v)n  &-m-\)(z)  “  •W-\>(z)] 


(22) 

(23) 


In  this  work,  we  always  take  m  a  positive  integer  or  zero.  For  later  con¬ 
venience,  we  also  record  the  formula 


VW  • (2‘> 

%. 


m,v 


Let  N  be  a  positive  integer  and  consider  that  solution  of  (21), 
call  it  cp„  /,(z)  with  ra  ^  N+2  such  that 

(25) 


(N)  _  (N)  _ 

^2, v  "  0  »  «PH+l,v  1  * 


Clearly  cp^N^(z)  is  a  linear  combination  of  the  solutions  (22)  subject  to 
the  conditions  (25)  and  we  readily  find  that 

~  z[Im+v^z^KN+2+v^Z^  +  6  KW-v^Z^IN+2+v^Z (26) 

in  view  of  the  Wronskian  relation 


Iv^Vl*2)  +  WZ>Kv(z)  =  i/z 


(27) 
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Suppose  that  we  are  given  the  normalization  relation 

03 

&(z)  =  H  ‘OfcVs/2)  • 

k=0 


(28) 


Put 


and  consider 


6 


(N),  \  V  fK ),  . 

(z)  =  Z_  <VPk,v(z) 


fc=0 


.(N).  . 

Vv(z) 


Q(z)tpl^u(z) 

9(N>(Z) 


m  s  N+l  . 


We  can  now  prove 


Theorem  5. 


lim 

N— >oo 


i^(z)  =  I  (z) 
mt-v  m+v  ' 


0  s  m  <  W+l  . 


(29) 


(30) 


(31) 


Proof:  Using  (26)  and  (29),  we  can  write 


.(N), 


e(z) 


= 


Wz)-(‘) 


N+m  IK+v+2^z^Km+v^Z'i 


%+v+2(z) 


N+3.  t  r z )  N+l 

k=0  %t-v+2'z'  k=0 


(32) 


and  the  result  follows  from  the  known  behavior  of  the  Bessel  functions 
for  large  order.  That  is. 


(  /2 )m+V 

■^hi+n/2)  ~  r(nt+v+l)  t+0(m  )]  > 

Vv(z)  =  |(z/2)"m'Vr(nt+v)[l+0(m-1)]  . 


(33) 

(34) 
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We  next  show  that  {pW(z)  ,  m  =  N+1,N,N-1,.  .  .  can  he  repre¬ 
sented  in  terms  of  a  generalized  hypergeometric  polynomial.  We  then  prove 
that  for  two  specific  choices  of  0(z)  ,  the  series  (29)  can  also  be  expressed 
in  terms  of  a  generalized  hypergeometric .polynomial;  and  further,  for  the 
two  choices  of  9(z),  respectively,  i(,  '(z)  and  the  rational  approxima¬ 
tions  tn(z)/hn(z)  and  en(z)/Sn(z)  »  respectively,  are  equal.  Actually, 
we  first  state  and  prove  theorems  for  the  Case  I  situation  in  some  detail. 

The  corresponding  theorems  for  Case  II  are  stated  and  proofs  are  omitted 
as  the  details  are  much  akin  to  their  Case  I  analogs. 

Another  choice  for  9(z)  previously  discussed  in  the  literature 
is  Iv(z)  •  We  call  this  Case  III  even  though  the  corresponding  ij^C2)  = 

^(z ) I^( z ) /<PqN^( 2 )  is  not  a  member  of  the  family  of  approximations  from 

which  Cases  I  and  II  were  derived.  We  defer  further  analysis  of  Case  III 
to  a  later  discussion  when  we  determine  closed  form  error  expressions  for 
all  the  cases. 

HYFERGECMETRIC  REPRESENTATION  FOR  cp^(z) 


etc. 
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(35) 


<p2n-2kte,v^z^  “  2—  (2m+2T|): 

m-0 


■e 


2(k+l-6  )(2n-k+l+v)~ 
z 


e6  -k,k+l+T]-e6  ,k-.2n-v-(l-6  )e,2n-k+2+v+e6  -6 

ini 


X  =  4/z^  ,  N  =  2n-6  ,  6  =  0  or  6  =  1  ,  e  =  0  or  e  =  1  , 


11  =  e6  +  (l-e)(l-6)  ,  (36) 

•where  6  =  0  or  1  according  as  N  is  even  or  odd,  respectively.  Also 

,  ,  .  .  ,2n-m+l-6  , 

„(»)(  Z)  =  (gA)  r(£n+2-6+v) 

m,v  r(m+v+l) 


r  m-l+61 

/  m-1+6 

m+6 

,1?  2  J 

3 

(  -n4“T~ 

,-n+— 

\  -2n-l+6 -v 

,m+v+l , -2n-l+6+m 

(37) 


where  [p]  is  the  largest  integer  s  p  . 


Proof:  By  induction:  The  Table  (35)  is  readily  developed  by  use  of  (21) 
and  the  starting  conditions  (25)  and  it  is  easily  verified  that  (36)  gives 
the  polynomials  listed  in  (35).  Put  (36)  with  6=1  and  e  =  0  in  (24) 
(in  (24)  replace  m  by  2n-2k  )  .  Then  after  some  algebra,  it  is  seen 
that  the  coefficients  of  like  powers  of  z  vanish,  which  proves  (36)  for 
6=1  and  e  =  0  .  To  get  (36)  when  6  =  $  =  1  ,  use  (21).  The  case  6=0 
is  similar  and  we  omit  the  details.  Finally,  (37)  is  just  a  special  case 
of  (36).  To  connect  these  two  equations,  we  set  2n-2k+e  =  m  and  choose 
€  =  0  or  1  according  as  m  is  even  or  odd,  respectively. 


Remark  1:  We  have  given  more  polynomials  in  (35)  than  are  necessary  for 
the  proof.  The  additional  entries  are  given  for  convenience. 


Remark  2:  If  in  (36),  e  =  6  =  1  ,  and  if  k 
and  v+1  ,  respectively,  then  we  get  (36)  with 
e  =  0  and  6=1,  and  if  v  is  replaced  by 
e  =  5  =  1  . 
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and  v  are  replaced  by  k*-l 
e  =  6  =  0  .  Again,  if  in  (36), 
v+1  ,  then  we  get  (36)  with 


11 


(41) 


Bk(»)  = 


(-)  (s-k):(vH), 


n-m  =  k  ,  n+m  =  s 


The  case  k  =  0  is  trivial.  Assume  k  >  0  .  Clearly  Bk(v)  is  a  poly¬ 
nomial  in  v  of  degree  k  which  vanishes  if  v  =  -u  ,  u  =  1,2,..., k  .  A 
straightforward  calculation  shows  that  br+bu  r  =  0  whence  h^(v)  also 
vanishes  if  u  =  l,2,...,k  .  Next  multiply  both  sides  of  (40)  by 
(k+-v+l)]t_^  .  Then  (40)  and  (41)  take  the  form 


hk(v)  =  Bk(v)  » 

,  =  v*  v*  v*  -  (-)r(2^v)r(^v)(-k)r(s^-n-v)rr(2^v) 
k  V  .  r  *  r  r.T(v+l)(-s)rr(k+l+v+r) 


*,  .  v-7 

Vv>  "  k!r(vn) 


Now  each  side  of  (42)  is  a  polynomial  in  v  of  degree  2k-l  and  b£(v) 
vanishes  for  u=  l,2,...,2k-l  .  Also  hg(v)  vanishes  for  u=  l,2,...,k  . 
If  v  =  -k-v  ,  v  =  l,2,...,k-l  ,  then  b£  =  0  for  r  =  0,1,..., v-1  and 
+  I’k-j  =  0  for  3  =  0,1, . . .  ,(k-v)/2  .  So  h£(v)  and  b£(v)  have 
the  same  zeros.  Further,  it  is  easy  to  see  that  the  coefficients  of 
on  both  sides  of  (42)  are  equal.  Hence  h£(v)  =  b£(v)  and  so  also 
hjc(v)  =  B)s(v)  for  all  values  of  v  which  proves  the  theorem. 


Theorem  8. 


*$(«>  =  (2/z) 


1-6  {-)n(v+l)n+l-6 


tn(z)  » 


N  =  2n-6  ,  6  =  0  or  6  =  1  . 


Proof:  By  induction:  Using  (4)  and  (37),  we  can  readily  verify  the  state¬ 
ment  for  n  =  0,  1  and  2  .  A  straight forward  analysis  shows  that  both 
sides  of  (44)  satisiy  the  same  recurrence  formula  which  is  easily  deduced 
from  (11). 
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From  Theorems  7  and  Q,  and  (32)  and  (37),  we  have 


Theorem  9. 


<»)  *n W 

^  .<*><*> 


(H)  (z/g)mtv  LI 

V',<z)  '  o(B)(z)  ’  r(»v+i) 


'  hn(z)  > 

|  _  m-1+6  |  /  m-1+6  m+6 
Ln - 2~\f  -n+  5  ,-n+-r- 


2*2 
-2n-l+6  -v,m+\r*-l,  -2n-l+6+m 


0F5(-2n-l+6-v;Z2/4) 


(45) 


Thus  the  result  produced  by  use  of  the  recurrence  formula  for  Iv(z) 
employed  in  the  backward  direction  together  with  the  normalization  relation 
(38)  and  the  rational  approximation  given  in  (2) -(5)  are  identical 


Case  II.  From  [l,  Vol.  2,  p.  45,  Eq.  (5)],  we  have  the  normalization 
relation 


n(z) 


=  (z/2]V  =  y-  (2k+2v)r(k+2y) 


r(v+i)  r(2v+i)k: 

k~  0 


v  i  0  > 


=  ez  =  I0(z)  +  2  Y_  Ik(z)  ,  v  =  0  . 
k=l 


(46) 


Here  to  avoid  confusion,  we  replace  0  by  n  in  the  notation  of  equations 
(28)  and  (29). 

Proofs  for  Theorems  10,  11  and  12  given  below  are  akin  to  those 
for  Theorems  7,  8,  and  9,  respectively,  and  we  skip  the  details. 


Theorem  10. 


nW(z) 


N+l 

z 


k=0 


(2k+2v)r(k+2v) 

r(2v+l)kl 
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(2^+2)^!  !  -N-1,N*2\H-3,1 


(IP-1).'  3  1  \v+3/2  ~1/2‘ 


(2v+2) 


IP1 


(nh-i):  giP-i^2) 


(47) 


Theorem  11. 

(N),  ^2s,+2Vl 

%v(z)  =  (N+i):  Vi(z)  » 

and  as  a  consequence  of  (44)  and  (48),  we  have 


<<48) 


(49) 


Theorem  12. 


(N)^\  =  _  9N+1^2) 

n^N\z)  WZ> 


K  '(*) 


m-1+61/  m-1+6  m+6 


b^r-V- 


«*— 5— ,-n*-  2 


/  /  ,m+v  z  2*3  \ 

^m,^2-1  (z/2)  e  _ \-2n-l+6  +m,nH-v+l,-2n-l+6  -m 

Vv^2)  =  q(n)(z)  =  r(m+v+l) 


1  1 


f  2n+l-6  /  -2n-3/2-v+6 


2z 


-4n-3-2v+26l 


(50) 


A  FURTHER  INTERPRETATION  OF  THE 
BACKWARD  RECURRENCE  PROCESS 


Fran  Watson's  treatise  on  Bessel  functions  [5,  p.  295],  we  can 


write 


I2n+2-6+v(z)  =  T2n-6  ,v+2(z)Iv^2]  "  T2n+l-6  .v+l^^v+l^2) 


(51) 


T2n+l-6  ,vU(z>  =  ‘PO,^2)  »  N  =  2n^  ,  6  =  0  or  6  =  1  ,  (S 

where  <PqN^(z)  is  defined  by  (37).  Consider  the  case  6=0  only.  We 
put  (51)  £n  the  form 


I\i(z) 


I2n+v<-2^z^ 


T2n+1, 


Further,  we  can  put 


Vfi 

Qn(z)  “  F  (z)  * 


/  ,2n+\»+l 

_  2>/2  T2n+l,v+l 

®n(z)  ~  r(2n+v*-2) 


.  .  ,2n+v+l 
_  T2n,v+2 

Fn(z)  "  r(2n+-v+2) 


Now  it  is  known  [5,  p.  302]  that  for  z  and  v  fixed,  z  /  0  , 


lim  En(z)  =  Iv(z)  ,  lim  Fn(z)  =  I^+^z) 


From  (44)  and  (52),  we  have 


En(z)  =  Ln(z)  >  Ln(^)  = 


(-)n(z/2)2nhn(z) 
n!  (n+v+2)n 


,  (56) 


where  \|rn(z)  and  h^z)  are  given  by  (3)  and  (5),  respectively,  with 
6=0.  But  for  z  and  v  fixed,  z  f  0  , 


lim  Ln(z)  =  1  , 


see  (9).  That  is, 


lim  tn(z)  =  lim  — — — r- — t— r  -  1 

r(w+l)h„(z) 

n— n->®  v  '  n'  '  n- 


(z/e)V  (z)  (z/2)Vcp^(z) 


r(v^l)0(N)(z) 


=  1M 
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A  similar  analysis  can  be  made  for  Fn(z)  .  Also,  a  like  study  can  be  done 
for  the  Case  II  scheme.  We  omit  the  details. 

Further,  the  above  shows  that  the  backward  recurrence  scheme  for 
the  computation  of  zlv(z)/lv+i(z)  is  the  same  as  the  well-known  truncated 
continued  fraction  representation  which  in  turn  is  the  same  as  the  main 
diagonal  Pade  approximation  for  this  function. 


ERROR  ANALYSES 


In  the  first  part  of  this  section  we  develop  closed  form  repre¬ 
sentations  of  the  error  in  i^^ ( z )  for  Cases  I-III  under  the  assumption 
of  exact  arithmetic.  This  type  of  error  arises  because  N  is  finite  and 
is  called  the  truncation  error.  From  each  analytical  representation  of 
the  error,  we  deduce  an  asymptotic  estimate  of  the  error  which  is  very 
realistic  and  easy  to  apply  in  practice.  The  results  for  Cases  I  and  II 
when  m  =  0  are  much  better  than  those  given  by  (10)  and  (19),  respectively. 
Further,  for  Cases  I  and  II,  if  z  ana  v  are  fixed  and  n  is  sufficiently 
large  with  respect  to  m  ,  the  relative  error  in  the  approximation  for 
I^^z)  is  essentially  independent  of  m  . 

An  analytical  formulation  of  the  round-off  error  is  developed  in 
the  latter  part  of  this  section  where  it  is  shown  that  this  source  of  error 
is  insignificant. 

We  now  turn  to  a  study  of  the  truncation  errors. 


Case  I .  Let 


where 


(N) 

Wz) 


«£!’<*>  -  W‘)  - 


is  given  in  (45). 


(59) 
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Theorem  13.  If  v  is  not  a  positive  integer  or  zero, 


0I^(-2n-l+6-v;z2/4)E^(z)  - 


(n+l):r(2n+2-6+v)  m+v 


,  (2)I  (z) 

f(2n+2-6+v)  sin  vrr  -m-vv  '  2n+2-8+v'  ' 


2(-)m+S(z/2)2nf2-6^ 

r(2n+2-6+v) 


I2n+2  -6  +  si  ( z  ) Kjjj+v  ( z ) 


.  .6  .  .  ,2n*-2-6+v 

(-)  tt(z/2)  _ 

+  f(2n+2-S+v)  sin  s>n  I2n+2-5+\/Z'Im+\FZ' 


Equation  (60)  can  be  rearranged  so  that  with  the  aid  of  L'Hospital's  theorem, 
we  can  get  a  representation  of  the  error  when  v  becomes  a  positive  integer 
or  zero.  We  do  not  give  this  result.  However,  for  arbitrary  v  ,  we  always 
have 


p*i(  2n  1+6  ,\)*z2/4)e1N)(z)  =  ^ ^Z/2^ - r_(n+l-6+y)  . 

(-2n-l  6  v,z  /4JEmjN;(zJ  (n+1) !r(2nf2-6+v)  Im+v^zaI2\n+2,-: 


n+6  -v 


2(-)m+6(z/2)2n+2'5+M  o((z/2)2n) 

+  r(2n+2-6+v)  I2n+2-6+y(z)Km+y(z)  +  r(2n*-2-6+\>)r(2iH-3-6+v)  ’ 

where  s  =  n-6+v(s=&>)  if  v  is  (is  not)  a  positive  integer  or  zero.  Clearly 
the  backward  recurrence  scheme  is  convergent.  Further,  for  n  sufficiently 
large,  n  »m  ,  the  relative  error  is  essentially  independent  of  m  .  For 
convenience  in  the  applications  we  record  the  formula 

(-)n(z/2)2rtf2r(ntl-6+v)  r  n  t 

=  ( ri+1)! r f2n+2 -6  +  v )  "  >] 


2(-)m+6(z/2)2nf2^+v 


r(2n+2-6+v)  I2n+2-6+\/z)Khi+v^ 
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ak+2n+2-6 -m  „4n+4-26  ■ 
2 


(-)m  5r(m+v)r(nH-\)-i-l) 


,2nV(2nt2-6+\))r(2n+3-6+\>) 


/  m  3-6  \  /  m  4-6  \ 

Vn  -  2  +  2  /  k  yn  ~  2  +  2  )k 

(2n+3-6+v)k(l-m-v)js.(2n+3-6-m)  k! 


,  k  i  0  .  (66) 


Thus 


(  -n+Sli^,.n+E± 

Vz)  =  2F3L  2  (  2  2 

\-2n-l+6  -'J,m*-s^l,-2n-l+6  +m 


,  (-)6nr(^^l)(Z/£)2rtfg~*  "m 


„ - W°>WW*>  •  «"> 


and  in  view  of  (23),  we  have  the  alternate  representations 


2n+2-6 -m 


Am( z )  ~ 


(-)6TTr(m+v+l)(z/2) 

r(2n+2-6+v)  sin  vrr 


Im+v^2)  [”  “  (~)  +  X2n+l-6+v^z 


2F3 


[n-lfi]  /_  n+S=±! ±,_n+5|L 


(-)  nr(nrt~v+l)(z/2) 
f(2n+2-6+v)  sin  vtt 


-2n-l+6  -v  ,m<- v+1,  -2n-l+6+m 
2n+2-5  -m 


I2n+2-6+v^z)[n  (-)"WZ)  +  Vv(zJ  * 

(68) 


As  a  remark  aside,  the  combination  (26)  and  (68)  yields  (37)  and  so  we  have 
an  alternative  proof  of  (37).  Using  (45),  (59)  and  (67),  we  can  write 
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oF^(-2n-l+6-v;z2/4)4^(z)  =  Vv(z)0F[n)(-2n-l+6-v;z2/4) 


.  (z/2)- 


m+\> 


r(m+v+l) 


^(nr+v+l ;  z2/4  )  qF1(  -2n-l+6  -v ;  z2/4 ) 


.  ,6  ,  w  .  . 2n+2-6  -m 

(-)  rrr(m+M*l)(z/2) 


r(2n+2-6+\j)  sin  \>n  -m-V  ;  2n+2-6+v 


]■ 


and  since 


0F£(-2n-l+6-v;z2/4)  -  ^(-Sn-l+S -v;z2/4) 


(-)nr(n+l-6+u) 

(n+l):r(2n+2-6+v) 


1F2 


1 

n+2 ,  -n+6 


(69) 

we  readily  find  the  first  part  of  (60) .  The  second  part  of  (60)  follows 
from  (23)  or  it  could  have  been  found  by  repeating  the  above  analysis  with 
the  second  equation  of  (68)  in  place  of  (67). 

Next  we  briefly  examine  the  situation  when  v  is  a  positive 
integer  or  zero.  With  v  =  r+g  ,  the  jFg  on  the  right-hand  side  of  (60) 
can  be  expressed  as 


F 
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(i 


iP-2,-n+5  -v 


1 

n+2 , -n+8 -v 


+ 


£ 


ism 


2k 


(n+2),  (-n+6 -v), 
k=n-6+r+l  'k 


=  Ji-6+r/l 

1  2  \n+2 , -nt6 


,n+l+6 


(-)  Tr(n+l) :  (z/2) 


2n+2-26+2r 


sin  nvr(2n+3-6+r)r(n+l-6+v)r(l-e) 


lF2 


I1 

\2n+3-6+r,l-e 


(70) 
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I 


The  first  term  on  the  right-hand  side  of  (70)  is  defined  when  e  =  0  . 

When  the  second  term  on  the  right-hand  side  of  (70)  is  multiplied  by  the 
coefficient  of  the  ^F2  in  (60)  and  the  result  is  combined  with  the  term 
involving  Im+v(z)l2n+2-6+v(z)  in  (60)>  ^  wil1  be  seen  ’bhat  one  can  pass 
to  the  limit  as  e— >0  .  The  final  expression  is  not  of  great  interest 
and  we  omit  further  details.  Equation  (61),  which  is  important  for  practical 
considerations  (see  later  numerical  example)  readily  follows  from  (60)  and 
the  above  remarks,  and  (62)  is  but  a  simplified  version  of  (61). 

Remark:  Let  v  ,  n  and  z  be  fixed  so  that  Ej^^(z)  is  a  function  of 
m  only.  Then  E^N^(z)  satisfies  the  recurrence  formula  for  cpm  ^(z),  see 
(21).  This  is  evident  from  (30)  and  confirmed  by  (60). 


Case  II.  Let 


Pw 

m,v 


(z) 


I 

m+v 


(z)  -  i^(z) 
m+v 


where  i£^(z) 


is  given  in  (50) . 


(71) 


Theorem  14.  If  neither  v  nor  v+5  is  a  positive  integer  or  zero,  then 


a"z  2n+l-6  ( -2n-3/2+6  -\>  \  (») 

'  1F1  \-4n-3+26-2v  2  / 


+  (^)gntg-6(v4)got^..-X^(z)  A  V>1 
+  (2n+2-6):(2n+2-6+2v)2rrf2_fi  2  2\2n+3-5 ,-2n-l+6-2v 


*) 


(-)\(z/2)dn+ 


-m-\/z  ^2n+2-6+v 


(2) 


r(2n+2-6+v)sin  vn 

which  is  the  same  as  the  right-hand  side  of  (60). 

If  is  a  positive  integer  or  zero,  call  it  r  ,  we  have 


(72) 
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where  t  =  v-^(t=2n+l-6+2\»)  if  v  is  half  a  positive  integer  (is  a  posi¬ 
tive  integer  or  zero)  and  where  t  =  ®  for  all  other  v  .  Clearly,  the 
■backward  recurrence  scheme  is  convergent. 

Further,  for  n  sufficiently  large,  n  »  m  ,  the  relative  error  is 
independent  of  m  .  For  convenience  in  the  applications,  we  record  the 
formula 

2n+2-6  ,  -z 

(N)  ~^2z^ _ ^  2^2n+2-6e  .  r  -i 

Fm,v(z)  =  (2n+2-6):(2n+2-6+2v)2n+2_6  Im+v(z) ^  +  °(n 


2(-)m+5(z/2)2n+2-6+v 

r(2n+2-6+v) 


I2iH-2-5+v^z)^n+v(z)  >  v^~2-(76) 


Proof :  Let 


J  z  /  m+\r*-^  I  \1  J"  _ z  /-2n-3/2+6 -v  I  _  V~  k 

"  |e  lFll2nH-2v+ll"27f  1®  lFl\-4n-3+26 -2v ‘2zV  I  f-  bkZ 


In  view  of  (1)  and  (12),  Bm(z)  =  A^z)  where  Am( z )  is  given  by  (64). 
Hence  b^  =  0  for  k  odd  and  b2^  =  ak  where  ak  is  given  by  (65)-(66). 
The  analysis  proceeds  as  for  Case  I  and  we  find  (72).  Notice  that  the 
right-hand  sides  of  (72)  and  (60)  are  identical. 


For  the  proof  of  (73),  let  Vm  v(z)  be  the  entire  first  term  on 
the  right-hand  side  of  (72),  that  is,  the  term  involving  the  gF2  • 
v  =  r+-^+e  ,  r  =  0,1,2,...  .  Then  we  can  write 


Wz)  =<Jz)(^l+e)2n+2 


2n+2-5+2r 

2k  , 

Z_  ckz  + 
k=r+l 


-e  [2^2- 2n+3-6,-: 


2n-l+6-2v! 


z 

k=2n+3-6+2r 


.2n+2-6  -z 

(2z>  e  Vv(z) 


^  Z;  ~  (2n+2-5 ):(2n+2-6+2v) 


2n+2-6 


(#-v) 

_ k _ 

*  Ck  (2n+3-6  )  (-2n-l+6-2v), 

K  K 


It  is  easy  to  pass  to  the  limit  when  e — *0  and  (73)  readily  follows. 
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Now  for  arbitrary  a  and  b  , 

x 

al  i(z)  +  bljjz)  =  (2/rrz)2(a  cosh  z  +  b  sinh  z)  , 

2  z 

Further,  from  (16),  with  n  replaced  by  2n+l-6  ,  see  also  (17),  and 
v  =  ,  gn(z)  reduces  to 

(n/2z)-ieZK2nt3/2_5(z)  . 


Use  these  data  with  r  =  -1  and  m  =  0  in  (73)  to  get  the  first  line  of 
(74).  Derivation  of  the  second  line  of  (74)  is  trivial  and  details  are 
omitted.  Since 

I_i(z)  =  (2nzfV(l+e-2Z)  , 

1  -z  (N) 

it  is  easy  to  show  that  (2rrz)2e  i  i/(z)  is  the  main  diagonal  Fade  approxima- 
-2z  ~2 

tion  to  1  +  e  .  When  allowance  is  made  for  a  change  of  notation, (74) 
is  a  previously  obtained  result  [l,  Vol.  2,  p.  74,  Eqs.  (34),  (35)]. 

When  v  =  r^e  ,  r  a  positive  integer  or  zero,  we  can  rearrange 
(72)  after  the  manner  of  the  discussion  surrounding  (70)  and  use  L’Hospital's 
theorem  to  get  the  limit  as  e  — >0  .  The  result  is  not  of  immediate  interest 
and  we  omit  details.  The  statements  (75)  and  (76)  are  readily  derived  and 
here  too  we  skip  details. 

(  n) 

fiemark:  F'  '( z )  with  n  ,  v  and  z  fixed  and  n  variable  satisfies  the 
-  m,v 

recurrence  formula  for  m  (z)  ,  see  (21).  This  feature  is  clearly  depicted 
by  (72).  In  both  (62)  and  (76),  the  term  involving  Kntfv(z)  is  of  lower 
order  than  the  term  involving  Im+N)(z)  .  Neglecting  the  former  term  in  each 
equation,  we  have 
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Theorem  15. 


Oz)  (-)n+1r(n+l^+v)(2nf2^)!22v-1(i)v>(Z/2)6eZ[L+0(n-1)] 
JN)7T=  (n*l):r(2n+2-6+2v) 


*<*>(0 

m,v' 


(-)n+1(i)v(2/2)6eZ 

- — -  [l+0(n-1)|  ,  v  i  4 


This  shows  that  there  is  little  difference  in  the  accuracy  of  the 
two  schemes  for  the  evaluation  of  Imfv(z)  .  Computation-wise ,  if  the 

backward  recursion  scheme  is  used,  Case  I  requires  less  operations  since 
the  associated  normalization  relation,  see  (38)  and  (39),  uses  the  sequence 

iaJN^(z)}  ,  k  =  0,2,4,...  ,  while  the  Case  II  normalization  relation,  see 
K,v  f  (N)  1 

(46)  and  (47),  employs  <^Cpk  ^(z)^  >  k  =  0,1,2,...  .  Also  to  get  I^(z) 

by  the  Case  II  scheme,  ez  must  be  evaluated.  On  the  other  hand,  if  |z| 
is  large,  r(z)  >  0  ,  one  often  wants  not  Iv(z)  ,  but  e"zIv(z)  .  The 
latter  is  automatically  furnished  by  the  Case  II  technique.  It  appears  that 
for  the  same  n  ,  the  Case  II  procedure  might  be  more  accurate  than  the 
Case  I  scheme  even  for  moderate  values  of  |z|  ,  R(z)  >  0  ,  in  view  of 
the  presence  of  ez  in  the  numerator  of  (78).  Also,  Case  II  is  favored 
when  R(v+6 )  <  0  .  Improved  information  cannot  be  derived  from  (78)  as 
the  estimate  is  for  fixed  m  ,  u  and  z  .  For  error  analyses  it  is  sug¬ 
gested  that  one  use  (62)  or  (76)  as  appropriate.  Further  discussion  is 
deferred  to  a  later  part  of  the  paper  where  numerical  examples  are  presented. 

If  z  is  pure  imaginary  and  v  is  real,  then  z-N,Iv(z)  is  real 
and  definitely  the  Case  I  procedure  is  better  than  the  case  II  scheme  since 
the  former  requires  real  arithmetic  while  the  latter  demands  complex  arith¬ 
metic  . 


If  only  Iv(z)  or  only  e_zIv(z)  is  required,  use  of  the 
rational  approximation  scheme  or  the  equivalent  backward  recursion  scheme 
demands  about  the  same  number  of  operations.  In  the  absence  of  a  priori 
estimates  of  the  error,  the  rational  approximation  scheme  employed  in  the 
following  fashion  is  preferred.  It  is  sufficient  to  consider  the  Case  I 
situation.  Compute  irn(z)  from  either  (3)  or  the  combination  (37),  (44), 
and  M*)  from  (5),  for  n  =  0,1  and  2  .  Compute  subsequent  values 


of  *n(z)  and  hn(z) 

tn(z)/hn(z)  ^  Vi 

one  requires  Ik*v(z) 


by  use  of  the  recursion  formula  (12).  Comparison  of 


(z)/hn^(z) 
or  e  *k+v 


affords  an  estimate  of  the  error.  If 
(z)  for  k=  0,1,2,. ,.,r  ,  then  obviously 


the  backward  recursion  scheme  is  highly  advantageous. 
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Case  III.  Consider 


G(N)(z)  =  I  (z)  -  i<N)(z)  , 
m,vv  '  nr*-v  m+\r  ' 

(M) 

W”  '  4tlM  ' 


■where  m'  '(z)  is  given  by  (37). 

m,\> 


Theorem  16. 


6+m  ,  (m-2) 

(N)  I2n+2-6+'j^Z^®0,v  ^ 

GW(Z)  =  -  - 

m,v  W.  . 

90,v(z) 

,  ,6+m  .  .  2n+2-m-6 

■  a  - WWV.W 


ra-l}/Jii  JSc£  \ 

pm(  2  ’  2  z2] 

2  3  \l-v-m,\>+l,l-m  / 

2  3  \-2n-l+6-v,  v+l,-2n-l+6 


,  0  <  m  s:  2n+l-6  ,  (81) 


and  it  is  clear  that  the  computational  scheme  is  convergent.  Further,  if 

we  treat  G'1  ( z )  with  n  ,  v  and  z  fixed  and  only  m  as  a  variable, 
m,v 

then  (-)mG^(z)  satisfies  the  recurrence  formula  for  m  (z)  ,  see  (21). 

m,v'  ^m,vv  '  v  ' 


(N)  ,  (N) 

Gm+2,v(z)  _  Gm,v(z)  _ 


n+v+l)  (N) 
z  Gm+l,v'z'  > 


6+1 

(N) ,  .  (N).  „  I2n+2-6+v'z ' 

Go>>  =  0  •  Gi>>  = - WT-, - 

«>0,v<z) 
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Finally,  for  convenience  in  the  applications,  we  record  the  formulas 


m+6  2n+4-25+2\>-m 

00,  ,  (-)  (g/2)  _ r(p^-v) 


m,\/z;  "  r(2n+2-6+v)r(2n+3-5+v)r(v+l) 


fm-l|/  m-1  m-2 

pItY  2  *  2 

3  \l-v-m,v+l,l-m 


[l+OCn-1)] 


/  \m+8  /  /oN2nt4-26+v-ri  , 
(N)  _  (-)  (2/2)  _ Tim 

r(2nK2-5+\>)r(2n+3-6+v) 


^  N.1  •  (84) 


Proof:  Using  (37),  (67)  and  (69),  we  find 


CPo!v(z)Gm!v(z)  =  '''knT1  I2i»^2-5+v<*)[Mz)l-m-v(*)-I-v(*)Vv(*)]  * 

and  in  view  of  (23) 

■  <->Tn+£!w2>  {*  [<->v>w>  -  • 


From  (26),  the  portion  in  curly  brackets  in  the  latter  equation  is 


(-)\(*~2)(z)  •  Hence  the  first  line  of  (81)  is  at  hand.  The  remainder 

of  (81)  follows  from  (37).  The  first  line  of  (81)  coupled  with  the  dis¬ 
cussion  surrounding  equations  (21)-(25)  produces  (82).  By  the  confluence 
principle,  see  [l,  V ol.  1,  p.  50l 


do,  ,  szm _ ui 

®o.v(z)  r(v+i) 


„  /  ,N+1  .  ,  f  N+ll  /  N+l  H 

2jz)  J-jlf— ,-g  ZS 

r^+1^  £  3  \-N-l-v,v+l,-N-l 


m-1 

(2/z)  r(N+v+2l 


0f1(v+1;z2/4)[i+o(it1)] 


N+'/*-l  r  i 

(2/z)  r(^v+2)iv(z)[i+o(r1)J  , 
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whence  (83)  and  (84)  readily  follow  with  the  aid  of  (33). 


Error-wise,  it  is  difficult  to  compare  Cases  I  or  II  with  Case  III 
without  some  simplifying  assumptions.  If  n»  m  ,  using  (33),  (62)  and 
(83),  we  get 

m-l|/  m-1  m-2 


*<*>(,) 

m,v'  ' 


(-)mfmf6(z/2)2mf2-26+2v-m(n+i):r(m^)  2Fs 


r ( 2n+3 -6  + V )r ( n+1-6  +V )r ( v+1 ) 
2n 


^ )  ^m+v 


-[l+0(n-1)] 


T(2n+3-6+v) 


.  6+1-v 
0(n  ) 


(86) 


and  so  Case  III  is  superior  to  Case  I.  Now  suppose  m  is  sufficiently 
large  so  that  in  (62),  the  second  term  dominates  the  first  term.  This  is 
certainly  the  case  if  m  =  2n+l-6-d  ,  d  «  n  ,  in  view  of  (33)  and  (34). 
Then 


a.) 


m,v' 


=  (2/z)6  [l+0(n'1)]  [l+OOn-1)) 


(87) 


and  under  these  conditions  there  is  little  to  choose  between  the  two  cases. 
Overall,  it  appears  that  Case  III  gives  better  accuracy  than  case  I.  How¬ 
ever,  for  Case  III,  one  must  know  Iv(z)  while  for  Case  I  no  such  knowledge 
is  required.  For  all  z  >  0  and  all  v  ,  0  s  v  s  1  ,  coefficients  are 
available  to  facilitate  the  rapid  evaluation  of  Jv(z)  and  Iv(z)  ,  see 
Cl,  Vol.  2;  6,7,8j.  (Actually,  much  more  information  is  given  in  these 
sources.)  All  of  this  can  often  make  the  Case  III  approach  rather  attrac¬ 
tive.  See  the  numerical  examples. 

Next  we  consider  the  round-off  error. 


Cases  I, II.  It  is  sufficient  to  trace  the  effect  of  a  given  round-off  error 
in  a  single  entry  of  the  table  generated  by  the  backward  recursion  process. 
Thus  let  Ctp  be  the  symbol  for  the  round-off  error  in  cp  .  Suppose  that 
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Thus  on  this  basis, 


Ai^(z)|s  a»|lv(z)  -  (z/2)V/r(v+l)|  (93) 


is  an  approximate  bound  for  the  round-off  in  a  single  entry  of  the  set  of 

numbers  generated  by  the  backward  recursion  process.  If  u>  is  the  maximum 

round-off  error  in  each  entry,  then  the  total  round-off  error  in  i'N'(z)  is 

m+\>  ' 

approximately  bounded  by  N  times  the  right-hand  side  of  (93).  Thus  the 
round-off  error  is  insignificant,  and  it  is  easy  to  estimate  the  number  of 
extra  decimals  which  must  be  carried  so  that  the  total  round-off  error  in 
the  process  lies  within  the  error  when  the  arithmetic  is  exact.  Equations 
analogous  to  (92)  and  (93)  for  Case  II  are  easily  derived  and  we  omit  details 


Case  III.  We  have 


,,00 


.  •  \  IX  J  f  \  _  ' 

Vj  z 


(94) 


Using  (26),  a  straightforward  computation  shows  that 

N,S.  .  .  . Y  2  f  .  6+Y+l  ,  .  1 

Am,v(z)  “  (")  z  [j2s+2-Y1'v(z'^2n+2-6+v'z'+(“'  -^2n+2-6+v'z'^2s+2-Y+v'z)J 


m+Y+1  (N)  (m-2) 

^2s+2-y,v^0,v  (2)>ra>0’ 


=  0  ,  m  =  0  . 


(95) 
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m+Y+1,  (N)  (N)  ,(m-2)  •> 

(N)  s  »(-)  Wz*P2s+2^v,v<z>Pn>  (z) 


= 


0,v  m,\>  (N)/  \ 

L  <>> 


Obviously,  this  is  nil  is  S  =  N  or  if  m  =  0  .  Now  using  (37),  we  find 


'N)  l-r\  i  ,2nt2s+3-6-V-m  .  „  , 

^23+2 -v,n/2^  _  (z/2) _  r(v-H)r(m<-v+l) 

{p(N)(z)tp(N)(z)  r(2n+2-6+v)r(2s+3-6+v) 


0,v  m,v 


[n.s.i!|zVj/  n+s+i^Y, .nts+i+^Y 


2n-l+6  -v,2s+3-v+^>-2ir^2s+l+6  -Y| 


F  P"^]  -n+f  z2\  fl 

\-2n-l+6 -v,v+l,-2n-l+6  j 


x  ,  m-1+6 1  /  m-I+6 

-t 

\-2n-l+6  -v  ,m+v+l, -2n-l+6+m 


and  it  is  clear  that  rounding  errors  are  insignificant.  Indeed,  if  all 
parameters  and  z  are  fixed,  then  the  round-off  error  tends  to  zero  as 
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FORMULAS  FOR  Jy(z) 


As  previously  remarked,  the  analyses  for  Iv(z)  hold  throughout 
the  cut  complex  z-plane,  -tt  <  arg  z  s  tt  >  and.  throughout  the  cut  complex 
v -plane  |arg  v|<  tt  >  although  it  is  sufficient  to  have  0  £  arg  z  s  tt/2  and 
R(v)  >  -1  .  Nonetheless,  we  indicate  how  to  get  results  for  Jv(z)  directly 
and  to  facilitate  use  of  our  findings,  it  is  convenient  to  restate  some  of  the 
key  equations.  We  omit  discussion  of  Case  II  since  it  requires  complex  arith¬ 
metic  to  generate  Jv(z)  which  is  real  when  z  and  \>  are  real.  In  any 
event,  the  reader  should  have  no  difficulty  in  establishing  the  Case  II  equa¬ 
tions  for  Jv(z)  once  it  is  observed  how  this  is  done  for  Case  I. 

All  developments  for  Jv(z)  are  readily  gotten  by  use  of  the 

equations 


-in/2 .  -i(m+v)n/2 

WZ®  }  =  6  Wz)  » 

(98) 

.  -in/2  ,  -i(m+s)n/2  (1) 

v,<«'  )  -  4"ie  Wz)  • 

(99) 

■  W‘) *  ‘W‘i  - 

(100) 

m  r  m  1 

=  (-)  esc  vn  [(-)  cos  Vtt  <Vv(z)-J_m_v(z)J  , 

(101) 

where  now  in  the  Jv(z)  analyses,  -tt/2  <  arg  z  s  3n/2  . 

It  is  convenient  to  introduce  the  following  notation.  Unless 
indicated  otherwise,  if  A  is  used  to  signify  some  function  or  equation 
in  the  developments  for  Iv(z)  ,  then  A*  is  used  to  signify  the  correspond¬ 
ing  function  or  equation  in  the  developments  for  Jv(z)  .  In  illustration 


^(z)  =  (z/2)  oFifv+lj-z2/*)  > 


(1)* 


and  both  J^^z)  and  Ym+V(z)  are  solutions  of  the  difference  equation 

Vv(z)  =  2(r  l}  Cl,v(z)  ■  *1+2, •  <21)* 
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Also 


/  /o  \  V 

<WZ)  =  ffm+v+TT  l+O^'1)  , 


.  .  ,-m-u 

Wz>  =  -  ~T  -  r^vjkm-1)] 


(34)* 


We  now  present  the  key  results  pertinent  to  Jv(z) 


Theorem  1*. 


(z/2)*  *n(z>  #  *  SJZ> 

Jv(z)  -  r(v+l)  *,  ,+  Bn(z)  »  ^(z)  ~  *  » 


V2) 


\(z) 


where  f*(z)  and  hn(z)  are  6iven  by  i^n(z)  and  hn(z)  ,  see  (3)-(5), 
provided  there  we  replace  z2  hy  -z2  ,  that  is,  replace  X  by  -X  . 
Further,  s£(z)  and  Sn(z)  are  both  given  by  (8)  and  r£(z)  =  (-J^RnC2)  > 
see  (10). 


Theorem  2.*  Both  t|fn(z)  and  hj(z)  satisfy  the  same  recurrence  formula 
(11)  if  there  we  replace  X  by  -X  . 

We  state  without  proof  the  following  equations. 


(N)*  _  (z/2)2n~m+1~5r(2n+2-5+v) 

'PmjV  2 '  T(m+v+l) 


X  2  o 


L_m-l+6 
L  2 


-n+m— +— ,  -n+^Q- 

2  2 

-2n-l-v+6  ,m+\H-l,-2n-l+6+m 


9  r(v+i) 


f-  (g^v)r(kH 

z—  r(v+i)k: 

k=0 


J2k*v(z) 


(38)* 


e(N)*(z)  =  (2/z)1"6  (v+1Vl-6  h*(z)  . 


(39)* 


33 


Theorem  9.* 


,  ,  (z)  **(z) 

,(N),  ,  * 

\  (z’  (N)*/  v  *,  \ 

9V  '  (z)  h  (z) 
n 


,  ,  VW*(*) 

.00,  , 

V,(; lz/2) 


m-1+61  /  m-1+6  , ,  m+6 

--^r\( 


»2n-l-v+6  ,m+\>+l,  -2n-l+6+m  1 


0F[n)(-2n-l+6-v;-z2/4) 


•  (45)* 


(N)*  (N) 

E  m,v(z)  =  JWZ)  -  Vv(z) 


in  J 

where  jm+v(z)  is  given  in  (45)*. 


Theorem  15.*  If  v  is  not  a  positive  integer  or  zero, 


y  2nt2  /,  ,  \ 


(  /  /ox2n+2-6+v 

ill  ^r(2n+2-6+v)  J-m-v(z)J2n+2-6+v(z) 

n(z/2)2n+2-5+V 

r(2n+2-6+v)  d2n+2-6+v'  '  nnVZ' 


.  ...  ,  ,2n+2-6+v 

(tt/2)(z/2) _ 

+  tan  vn  r(2n+2-6+v)  Jm+v(z^J2n+2-6+v(z/ 


(60)* 


Equation  (60)*  can  be  rearranged  so  that  with  the  aid  of  L'Hospital's 
theorem,  we  can  get  a  representation  of  the  error  when  v  becomes  a  posi¬ 
tive  integer  or  zero.  This  result  is  omitted.  However,  for  arbitrary  v  , 
we  always  have 
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r 


ftieorem  16.* 
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NUMERICAL  EXAMPLES 


Let 


N=5,n=3,6  =  l,z=2/3,\>=l/3. 

Values  of  ,  9^(z)  and  fi^(z)  are  given  in  the  table  below. 

(N) 

m 

6  1 

5  19 

4  305 

3  3984 

2  40145 

1  2  84999 

0  11  80141 


Since 

(z/2)V  -2/3 

Vi'S  =  0.77645  82114,  e  '  =  0.51341  71190  , 

r(v+i) 

the  Case  I  and  Case  II  approximations  are 


0.84272  08930  and  0.84272  10326 
respectively.  To  10  decimals 


Il/3(2/3)  =  0.842772  08819  . 

Thus  the  errors  in  the  Case  I  and  Case  II  approximations  are 

-7  -6 

-0.111-10  and  -0.151*10  , 

respectively.  Using  (62)  and  (76)  each  with  0(n_1)  and  the  term  involving 
Km+^(z)  neglected,  the  roproximate  Case  I  and  Case  II  errors  are 


0^(z)  =  880  75120/81 
n^(z)  =  1  38952  97360/6561 
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Here  the  entries  in  the  Im(z)  column  are  correct  for  the  number  of 
decimals  given. 

Using  the  first  lines  of  (62)  and  (76),  each  with  0(n_1)  and 
the  term  involving  Kjn+^z)  neglected,  the  approximate  relative  error 
for  Cases  I  and  II,  respectively,  are  -0.U6.10"3  and  -0.537-10"4,  respec¬ 
tively. 


In  the  table  below,  we  record  the  approximate  errors  obtained  by 
use  of  (62)  with  0(n-1)  omitted  for  m  =  6  and  5  and  by  use  of  (21),  see 
the  remark  following  Theorem  13,  for  tie  lower  values  of  m  .  This  is 
called  Case  I,  (62)-(21)  in  the  table.  We  also  present  the  analogous  Case  II, 
(76)-(21)  data.  In  each  instance  known  tabular  values  of  1^(2)  and  lj(2) 
were  used.  In  practice,  we  suggest  using  (34)  or  the  lead  term  of  the 
uniform  asymptotic  expansion  of  K^+<J(z)  developed  by  Olver  [9] .  For 

l2m2-6+\)(z)  >  use  (33)  or  the  lead  term  in  the  uniform  asymptotic  expansion 

for  this  function  which  is  also  given  in  the  source  just  cited.  We  also 
suggest  that  computation  of  the  gamma  functions  be  simplified  as  follows. 

With  |R(or)|  <  1  and  r  a  positive  integer,  we  have 


r(r+»+l)  =  r! 


r(rtar»-lj  _ 
r(m)  - 


rlr*  [l+0(r-1)] 


and  for  r  sufficiently  large,  we  neglect  0(r-1)  .  The  approximation  is 
of  course  superfluous  if  y  =  0  .  If  9  =  ±  j  ,  the  approximation  may  still 
be  used  though  known  tables  of  the  gamma  function  for  half  an  odd  integer 
may  be  preferred  [lO ] .  If  more  precise  values  of  the  gamma  functions  are 
required,  see  [lli  . 


Approximate 

Error 

m 

Case  I,  (62)-(21) 

Case  II,  (76)-(21) 

0 

-0.264-10-3 

-0.114-10"3 

1 

-0.184.10-3 

-0.765-10'4 

2 

-0.797 -lO'4 

-0.373.10"4 

3 

-0.242.10-4 

-0.110-10"4 

4 

-0.723- lO-5 

-0.410- 10"5 

5 

0.475.10-5 

0.536-10"5 

6 

-O.3iO.lO-4 

-0.309.10"4 

For  a  final  example,  we  illustrate  Case  III  using  the  data  of  our 
second  example.  We  get  the  following  numbers. 
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m 

i <"><»> 
m 

tW(.) 

m+v 

Approximate 

Error 

& 

00 

1 

CO 

(84) 

0 

2.27858 

5302 

0 

0 

0 

1 

1.59063 

6693 

0.162*10“® 

0.121*10“6 

0.276*10“® 

2 

0.68894 

8609 

-0,162 *10"® 

-0.121*10”® 

-0.276*10“® 

3 

0.21273 

9475 

0.484*10“® 

0.363*10“® 

0.551*10“® 

4 

0.50730 

1826 -10'1 

-0.161*10“® 

-0.121*10“®  -0.165*10“® 

5 

0.98187 

4502 *10“2 

0.693* 10“5 

0.520*10“® 

0.661*10“® 

6 

0.16364 

5750* 10“2 

-0.363* 10"4 

-0.272 *10“4 

■0.331.10-4 

In  the  above,  the  approximate  error  (82)- (83)  means  that  G^+i(z)  is 

approximated  by  (83)  with  m  =  1  and  0(n“^)  neglected,  and  subsequent 
approximate  values  of  the  error  are  found  by  use  of  the  recursion  formula 
in  (82),  Use  of  the  recurrence  formula  in  this  fashion  is  stable  as  the 
magnitude  of  the  error  in  an  increasing  function  in  m  .  Also  Eq.  (84) 
means  this  equation  with  0(m"-*-)  and  0(n"^)  neglected. 


A  measure  of  the  accuracy  of  the  three  schemes  treated  can  be  had 
by  use  of  normalization  relations.  Thus  if  the  Case  III  procedure  is  employed, 

(N) 

then  (38)  and  (46)  with  lj£+v(z)  replaced  by  ik+u(z)  are  available  as 

checks.  Similarly,  equations  (46)  and  (38)  are  available  as  checks  for  the 
Case  I  and  Case  II  techniques,  respectively.  For  some  other  useful  normali¬ 
zation  relations,  see jjL,  Vol.  2,  pp.  45,  4g] . 

Analyses  of  the  error  in  the  backward  recursion  process  for  the 
solution  of  a  general  second  and  higher  order  linear  difference  equation 
have  been  given  by  a  number  of  authors.  Some  authors  have  studied  the 
case  of  Bessel  functions  directly.  We  make  no  attempt  to  survey  the  various 
contributions  here.  Pertinent  references  are  given  by  Wimp  [4l  .  Suffice  it 
to  say,  none  of  the  analyses  have  the  precision  and  simplicity  of  those 
developed  in  the  present  paper.  We  deliberately  chose  N  and  as  a  conse¬ 
quence  n  small  (N=5,n=3 )  in  our  numerical  examples  to  put  our  asymptotic 
estimates  under  a  severe  test.  The  efficiency  and  realism  of  our  error 
formulas  is  manifest. 


CONCLUDING  REMARKS 


It  appears  that  the  techniques  developed  here  for  the  Bessel  func¬ 
tion  Iv(z)  can  be  extended  to  analyze  more  general  second  and  higher  order 
difference  equations.  In  particular,  it  would  be  useful  to  have  analogous 
results  for  gF;|_(a,b;c ;z )  and  its  confluent  forms.  This  we  intend  to  do 
in  future  papers. 
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